Zero-bias anomalies and boson-assisted tunneling through quantum dots 
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We study resonant tunneling through a quantum dot with 
one degenerate level in the presence of a strong Coulomb re- 
pulsion and a bosonic environment. Using a real-time ap- 
proach we calculate the spectral density and the nonlinear 
current within a conserving approximation. The spectral den- 
sity shows a multiplet of Kondo peaks split by the transport 
voltage and boson frequencies. As a consequence we find a 
zero-bias anomaly in the differential conductance which can 
show a local maximum or minimum depending on the level 
position. The results are compared with recent experiments. 
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Transport phenomena through discrete energy levels 
in quantum dots have been studied by perturbation the- 
ory [0J^] and beyond H-§|. In general, resonant tunnel- 
ing phenomena and Kondo effects in nonequilibrium be- 
come important, which have been measured recently by 
Ralph & Buhrman ||. In metallic islands, the Coulomb 
blockade is strongly influenced by inelastic interactions 
with bosonic degrees of freedom, such as fluctuations 
of the electrodynamic environment or applied time- 
dependent fields [||. The study of inelastic interactions 
in quantum dots with few levels has started only recently, 
either for the nondegenerate case ||Jl(| or more general, 
in the presence of time-dependent fields and Coulomb 
blockade (^,|^. In earlier work we have studied the in- 
fluence of bosonic fields in the noncquilibrium Anderson 
model in the perturbative regime pd[ | and found resonant 
side peaks in the Coulomb oscillations. 

The purpose of the present letter is to investigate the 
influence of external quantum-mechanical fields on trans- 
port phenomena through ultrasmall quantum dots at low 
temperatures and frequencies (compared to the intrin- 
sic broadening of the resonant state in the dot). This 
requires a description of the Kondo effect, generalized 
to nonequilibrium situations and including coupling to 
bosonic fields. For the nonperturbative treatment of the 
tunneling we apply a real-time, noncquilibrium many- 
body approach developed recently Jl2|,[l3| to a quantum 
dot with one level and spin degeneracy M. For M > 2 
and low lying dot level e we obtain the usual Kondo peaks 
at the Fermi levels /j, a of the reservoirs jjj. However, the 
emission of bosons causes additional Kondo singularities, 
for a one mode field at /i Q + najs {n = ±1, ±2, . . .). 

Furthermore, we will analyze the effect of the singu- 
larities in the spectral density on the differential conduc- 



tance as function of the bias voltage. For a low lying 
level we obtain the well-known zero bias maximum [Q-^) , 
whereas for a level close to the chemical potentials of the 
reservoirs we find a zero bias minimum. The coupling 
to bosons gives rise to satellite anomalies, which can be 
traced back to the corresponding satellite peaks in the 
spectral density. In a certain range of gate voltages, for 
M = 2 and in the absence of bosons, we find that the 
temperature and bias voltage dependence of the conduc- 
tance coincides with recent measurements of zero-bias 
minima in point-contacts ff4|| . Therefore, in addition to 
Refs. [p^|-[r^| , we propose here another possible interpre- 
tation of this experiment. 

We consider a dot containing only one energy level 
with degeneracy M connected via high tunnel barri- 
ers to reservoirs of noninteracting electrons. We, fur- 
thermore, include a coupling to bosonic modes repre- 
senting phonons, photons or fluctuations of the electro- 
dynamic environment. Our model Hamiltonian reads 
H = H + Ht, where H describes the decoupled sys- 
tem and Ht the tunneling between leads and dot. We 
write H = H R + H D where H R = J^kaa e kaal aa a k aa 
refers to the reservoirs (<r and a are spin and reservoir 
indices). Furthermore, (h — k R = 1) 
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H D = e N + U Y,<r<r' 

+ J2 q u q d\d q + NJ2 q g q (d q + d\) 



(1) 



describes the isolated dot with M spin degenerate levels 
at position eo, Coulomb repulsion [7 , bosonic modes u q 
and electron-boson coupling g q . The number of particles 
on the dot with spin a is denoted by n a = c^c a , and 



N 



Finally, the tunneling term is given by 



h.c). 



This Hamiltonian can be rewritten after a unitary 
transformation [Q defined by V = exp(—iNip) and (p = 
*E,GfeM)(4 " d i)- We g et H = VHV- 1 = H Q + H T , 
where H = H R + H D , H D = eN + J/J) , n a n a > + 
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rjf+h.c). Due to 



the electron-boson interaction the level position and the 
Coulomb repulsion are renormalized, e = eo — Egffg/^g 
and U = Uo — ^'Yli q 9 q l UJ q^ an d t ne tunneling part con- 
tains now phase factors e ±Iv . 

In lowest order perturbation theory the rates for tun- 
neling in and out of the dot to reservoir a are 



J dE'j^(E')P ± (E — E') 



(2) 
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where 7±(J5) = 1 /( 27r Wc t (E)f±{E) is the classical rate 
without bosons, T a (E) = 27r£ fc \ T k\ 2§ ( E ~ e k a ), and 
fa (E) is the Fermi distribution of reservoir a with chem- 
ical potential fi a while /~ (E) = l — f+(E). Furthermore, 
P^ (E) describes the probability for an electron to absorb 
(for P + ) or emit (for P~) the boson energy E. It is § 

P±(E) = -!- / dte tEt < j<e(°) e -M±t) >q ^ 
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where < ... >o denotes the expectation value with re- 
spect to the free boson Hamiltonian. The classical rates 
together with a master equation are sufficient in the per- 
turbative regime V = J2 a ^a < ^ @- In this letter 
we are interested in temperatures and frequencies which 
are of the order or smaller than the intrinsic level broad- 
ening r, which requires a nonperturbative treatment in 
r. To achieve this we use a real-time technique devel- 
oped in Jl^,[lj| which provides a natural generalization 
of the classical and cotunneling theory to the physics of 
resonant tunneling. For details we refer to these papers. 
Here we only sketch the derivation and quote the results. 

We develop a diagrammatic approach by expanding in 
the tunneling Hamiltonian Ht- Since Hq contains in- 
teraction terms, this can not be done by usual Green's 
function techniques since Wick's theorem does not apply. 
However, we can use it with respect to the field operators 
of the reservoirs, since Ho is bilinear in these operators. 
As an example let us consider the reduced density op- 
erator of the dot. We assume that the reservoirs and 
the boson bath remain in thermal equilibrium. On the 
other hand, we want to study the nonequilibrium time 
evolution of the dot. An effective description in terms of 
the dot degrees of freedom can be derived by expanding 
all propagators in Ht and tracing out the reservoirs by 
applying Wick's theorem for them. A matrix element of 
the reduced density operator can then be visualized as 
shown in Fig. (Q) . The forward and the backward propa- 
gator (Keldysh contour) are coupled by "tunneling lines" 
associated with the junctions to each reservoir a. Each 
tunneling line with energy E represents the rate ^(E) 
if the line is directed backward with respect to the closed 
time path and 7" (E) if it is directed forward. Because of 
Fermi-Dirac statistics, we get a factor —1 if two tunneling 
lines cross each other. The tunneling lines are associated 
with changes of the state of the dot, as indicated on the 
closed time path. For large Coulomb repulsion U we can 
restrict ourselves to states with N = 0, 1. The coupling 
to the bosons is introduced by connecting all vertices in 
all possible ways by boson lines with a certain energy E 
(the direction can be chosen in an arbitrary way). The 
rule for the contribution of the boson lines is exactly the 
same as for the reservoir lines except that we have to 
replace 7^ by P^. Finally, as in the case for metallic is- 
lands jl2| , we have to associate to each tunneling vertex 
at time t on the contour a factor exp i AE t where AE is 
the difference of the out- and incoming energies. If the 



vertex lies on the backward propagator it acquires a fac- 
tor — 1. Analogous graphical rules hold for the Green's 
functions of the dot as well, the only difference is the 
occurrence of external vertices. 

In leading order, we include only boson lines between 
vertices which are already connected by tunneling lines. 
This simply amounts to a dressing of the tunneling lines 
7^7, and the diagrams look identical to those without 
bosons. The approximation, while neglecting many dia- 
grams, describes well the spectral density of the dot at 
resonance points. The reason is that position and value 
of the peaks of the spectral density are determined by a 
self-energy a (see Eq. (^)) which is calculated here in low- 
est order perturbation theory in F including the bosons. 
Higher orders are small for high tunnel barriers. 

Similar to the case of metallic islands Jl2|,[l3| we pro- 
ceed in a conserving approximation, which takes into ac- 
count non-diagonal matrix elements of the total density 
matrix up to the difference of one electron-hole pair ex- 
citation in the reservoirs. The difference in the case of 
a single level quantum dot is that we have now M pos- 
sibilities for the occupied state. The analytic resumma- 
tion of the corresponding diagrams yields for the tran- 
sitions between N = 0(1) and 1(0) the rates = 
A / d^7 ± (^)|i?(£;)| 2 . Here R{E) = [E - e - a-(E)}- 1 
defines a resolvent with broadening and energy renor- 
malization given by the self-energy 



a(E) 



dE 



, M 7 +(E')+^-(E') 
E-E' + i0+ 



(4) 



where 7 ± = £ Q it and A" 1 = JdE\R{E)\ 2 . 

In the classical limit T << T we recover for X^ the 
classical rates r y ± . The stationary probabilities Pq and 
Pi for an unoccupied or occupied dot state follow from 
the kinetic equation which uses the rates as input P\T,~ — 
P ME+ = 0. Together with P + P x = 1 we obtain P = 
J dEj- (E)\R(E)\ 2 and P x = M J dEj+ {E)\R{E)\ 2 . 

Summing equivalent diagrams for the real-time Green's 
functions of the dot we obtain the spectral density p = 
{G<-G>)/{2m) 

p(E) = fdE'J2 Y(E')P- r (E> E)\R(E')\ 2 (5) 

J r=± 

and for the current I a flowing into reservoir a 

I a = e2nMj2 [ dE' ' ]T r 7 ~' \E')YAE')\R{E')\ 2 . (6) 

a' J r=± 

For the special case of two reservoirs a = L/R and 
constant level broadening r = Tl — T/j the current 
I = II = —Ir can be written as 



M f 

I = e-T J dEp(E)[f+(E)-f+(E)] 



(7) 
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Our results satisfy all sum rules together with current 
conservation, and one can prove particle-hole symmetry 
in the case M = 1. 

The difference to other approaches in the M = 1 case 
P,[To| is clearly displayed by the effect of the self-energy 
a(E) which determines via the resolvent R(E) the posi- 
tion of the maxima of the spectral density (||) . In all pre- 
vious works, a(E) has been approximated by a constant. 
We find that the energy dependence of <j{E) cannot been 
neglected if the temperature T and the typical frequency 
u>b oi the bosons are smaller than T. To derive this an- 
alytically we consider from now on a one-mode environ- 
ment (Einstein model) with boson frequency w q = lub- 
Defining g = E q 9 2 q /^B we obtain P ± {E) = £„p„<5(£± 
nuj B ), where p n = e -^ 1+2N ^ e nuJB / 2TB I n {2gN Q e UB / 2TB ) 
is the probability for the emission of n bosons with fre- 
quency ll>b- Here, No is the Bose function, /„ denotes the 
modified Bessel function. The temperature of the boson 
bath is Tb- In real experiments it can be different from 
the electron temperature T. Using ^) and <Q) we obtain 

Rea(E) = £(Mp„ - P -n)^ [in (||) 

and Ima(E) = -TrJ2 n Pn[Mj+(E + nuj B ) + J~(E - 
nus)]. Here <3/ denotes the digamma function, and we 
have chosen in the energy integrals a Lorentzian cut-off 
at Ec- The real part of a depends logarithmically on 
energy, temperature, voltage and frequency. These log- 
arithmic terms are typical for the occurrence of Kondo 
peaks and do not cancel for M > 2 or p n ^ P- n - Hence 
we anticipate logarithmic singularities not only for the 
degenerate case but also for a single dot level without 
spin since the probabilities for absorption and emission 
of bosons are different. This is an important difference to 
the case of classical time-dependent fields [|] where both 
probabilities are equal. 

Fig. (^) shows a typical series of pictures for the spec- 
tral density at different voltages for a low lying level e. 
Without bias and M — 2, we obtain the usual Kondo 
peak near the Fermi level (which we choose as zero en- 
ergy). Due to emission of bosons, there are now addi- 
tional resonances at multiples of lub ■ For finite bias volt- 
age, all peaks split and decrease in magnitude. 

The resonances in the spectral density can be seen most 
pronounced in the nonlinear differential conductance as 
function of the bias voltage V. Fig. (||) shows the dif- 
ferential conductance for a low lying level e. As usual 
we find a zero-bias maximum [|]-[6| since the splitting of 
the Kondo peak leads to an overall decrease of the spec- 
tral density in the energy range \E\ < eV (see inset of 
Fig. (|)). Due to emission of bosons we observe also a set 
of symmetric satellite maxima. They can be traced back 



to the fact that pairs of Kondo peaks can merge if the bias 
voltage is given by multiples of the boson frequency (see 
Fig. This gives rise to pronounced Kondo peaks 

at E = ±eV/2 and thus to an increase of the spectral 
density with bias voltage near these points. 

The differential conductance for e near zero energy is 
shown in Fig. (|]) with and without bosons. Surprisingly 
we find that the whole structure is inverted compared to 
the e < case and we find a zero-bias anomaly although 
the Kondo peak at zero energy is absent. The contribu- 
tions of sequential and cotunneling lead only to an overall 
shift of the differential conductance without any interest- 
ing structure. This shows clearly that the influence of the 
logarithmic terms in a{E) is still important. They lead 
to an overall increase of the spectral density near zero 
energy with bias voltage. In the presence of bosons we 
obtain satellite steps at \eV\ — itiujb- 

The occurrence of zero-bias minima is well known for 
Kondo scattering from magnetic impurities ||l9[| . Here we 
have shown that zero-bias minima can also occur by res- 
onant tunneling via local impurities if the level position 
is high enough to enter the mixed valence regime. We 
have also compared the scaling behavior of the conduc- 
tance as function of temperature and bias voltage with 
recent experiments of Ralph & Buhrman (m) (see insets 
of Fig. ([|)). The coincidence is quite remarkable. The 
explanation of this experiment, either interpreting it as 2- 
channcl Kondo scattering from atomic tunneling systems 
]f5| , [T6[ or by tunneling into a disordered metal [jlTj , is still 
controversial. The mechanism described in this work of- 
fers another possibility although the magnetic field de- 
pendence of the experiments remains unexplained. 

Finally, we also investigate the differential conductance 
at fixed bias voltage as function of the position of the dot 
level, which experimentally can be varied by changing 
the gate voltage coupled capacitively to the dot. Fig. (g) 
shows the classically expected pair of peaks at |e| = eV/2 
together with satellites between the main peaks (due to 
emission and absorption) and peaks for |e| > eV/2 (only 
due to absorption of bosons). The imaginary part of <r(E) 
gives rise to a classically unexpected asymmetry of the 
peak heights. The peak at e = eV/2 is larger than the one 
at e = -eV/2 since \Ima(E)\ = n\M-f + (E) + J~{E)\ is 
always smaller for higher energies (except for the M = 1 
case where particle-hole symmetry holds). This demon- 
strates a significant effect due to the broadening of the 
spectral density by quantum fluctuations. 

In conclusion, we have studied for the first time low- 
temperature transport in the noncquilibrium Anderson 
model with bosonic interactions. A one-mode environ- 
ment yields new Kondo resonances in the spectral density 
which can be probed by the measurement of the nonlinear 
differential conductance. We have shown that both the 
gate and bias voltage dependence is important. Quan- 
tum fluctuations due to resonant tunneling yield zero- 
bias anomalies as function of the bias voltage, which can 
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be changed from maxima to minima by varying the gate 
voltage. We found similarities to recent experiments. 
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FIG. 1. A diagram showing various tunneling processes: 
sequential tunneling in the left and right junctions, a term 
preserving the norm, a cotunneling process, and resonant tun- 
neling. 
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FIG. 2. The spectral density for T = T B = O.Oir, e = -4T, 
g = 0.2, lob = 0.5r and E c = 100r at different voltages. For 
V — there are resonances at multiples of u>b, which split for 
finite bias voltage. 

FIG. 3. The differential conductance vs. bias voltage for 
T — Tb — O.Oir, e = -4r, cub = 0.5r and E c = 100P. The 
curves show a maximum at zero bias and satellite maxima 
at multiples of lob for a finite electron-boson coupling. Inset 
(g = 0) : increasing voltage leads to an overall decrease of the 
spectral density in the range \E\ < eV, which explains the 
zero-bias maximum. 

FIG. 4. The differential conductance vs. bias voltage for 
T = T B = O.Oir, e = 0, lob = 0.5V and E c = 100F The 
curves show a minimum at zero bias and steps at multiples 
of for a finite electron-boson coupling. Left inset: the 
rescaled curves for g = at different temperatures collapse 
onto one curve. Right inset: The temperature dependence of 
the linear conductance (solid line) coincides with experimen- 
tal data from Jl4|] (triangles). 

FIG. 5. The differential conductance as a function of e for 
T = 0.25F eV = 30r, g = 0.3, lo b = 5T and E c = 500F 



4 



